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摘要 : 为 了 确定 栖息 地 的 破碎 化 和 片断 化 引起 中 国 大 饮 的 地 理 分 化 和 遗传 结构 变异 ， 本 文 测定 了 来 自 广 
西 、 河 南 、 陕 西 和 湖南 4 个 地 理 种 群 的 28 条 大 鲍 的 mtDNA D-loop 基因 全 序列 。 根 据 这 4 个 地 理 种 群 的 地 理 分 
布 , 分 成 珠江 单元 (广西 种 群 )、 黄 河 单 元 (河南 种 群 ) 和 长 江 单 元 (湖南 和 陕西 种 群 )。 通 过 Clustal X. 
MEGA2.0、DnaSP4.0、Arlequinl1. 1 分 析 发 现 ， 全 序列 长 度 为 771 bp， 其 中 64 个 多 态 性 核 苷 酸 变异 位 点 ， 占 全 
部 碱 基数 的 8.2690. TEHSCRURIUA A 37 6781 2 个 ， 插 人 /缺失 11 个 。27 个 单 倍 型 间 的 序列 差异 平均 为 1.32%, 
3 个 单元 的 单 倍 型 多 样 性 指数 和 核 苷 酸 多 样 性 指数 值 都 偏 低 ， 而 且 珠江 单元 的 这 两 个 指数 值 都 低 于 长 江 和 黄河 
两 个 单元 。 珠 江 单 元 和 黄河 、 长 江 单元 之 间 分 化 程度 显著 (P < 0.001)， 而 长 江 和 黄河 单元 之 间 差 异 不 显著 
(P >0.05)。 地 理 单元 内 分 化 程度 占 99.31% ， 而 单元 间 只 有 0.69%， 表 明 遗 传 差异 主要 发 生 在 单元 内 ,而 且 
各 地 理 单元 之 间 的 基因 流 较 频繁 。 构 建 的 NJ 树 和 MP 树 显示 ，27 个 单 倍 型 呈现 一 种 混杂 的 分 布 格局 ， 并 未 分 
成 代表 3 个 地 理 单元 的 聚合 群 。 
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Genetic Structure and Geographic Subdivision 
of Four Populations of the Chinese Giant 
Salamander (Andrias davidianus ) 
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Abstract: In order to identify the relationship between geographic subdivision and population genetic structure due to 
habitat loss and fragmentation, the complete sequences of the mitochondrial control region from 28 Chinese giant salaman- 
ders collected from Yangtze River Region (Hunan and Shaanxi), Yellow River Region (Henan) and Pearl River Region 
(Guangxi) were sequenced. According to the analysis of Clustal X, MEGA2.0, DnaSP4.0 and Arlequinl.1, 64 nu- 
cleotide sites were variable, including 6 transitions, 2 transversions and 11 insertions/deletions in the complete D-loop 
gene sequence (771 bp). A total of 27 haplotypes were identified and the average genetic distance among all these haplo- 
types was 1.32% . Furthermore, the results showed that all three geographic regions exhibited low levels of nucleotide di- 
versity and haplotype diversity, in which the two indexes of Yangtze River were the lowest. The results also presented that 
there was significantly distinctive geographic phylogeography between Pearl River Region and Yangtze River Region, be- 
tween Pearl River Region and Yellow River Region (all P « 0.001), but no obvious subdivision between Yangtze River 
Region and Yellow River Region (P » 0.05). In addition, and the genetic variation within and between regions was 
99.31% and 0.69% respectively, which means the genetic variation mainly happens within the geographic region and 
gene flow among the geographic regions was comparatively frequent. At last the phylogenetic tree analysis (NJ method and 
MP method) demonstrated that all the haplotypes could not be clustered to three geographic clades, but appear disordered 
distribution pattern. 
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兰 椎 动物 栖息 地 的 破坏 和 过 度 利 用 使 许多 物种 
的 种 群 数量 极度 下 降 ， 并 影响 其 生存 (Birstein, 
1993; Doukakis et al, 1999; Vos & Chardon, 
1998) 。 两 栖 类 的 扩散 能 力 比较 差 且 一 般 都 有 归 家 
行为 ， 故 栖息 地 的 恶化 、 种 群 大 小 剧烈 变动 更 易 导 
致 其 种 群 遗 传 多 样 性 的 丧失 ， 而 遗传 漂 变 、 近 亲 繁 
殖 的 不 断 增加 也 进一步 加 剧 了 两 栖 类 动物 的 濒危 程 
度 (Blaustein et al，1994; Seppa et al, 1999), B 
此 ， 栖 息 地 破碎 化 对 两 栖 类 遗传 多 样 性 和 地 理 种 群 
差异 的 影响 成 为 研究 的 热点 (Reh & Seitz, 1990; 
Hitchings & Beebee, 1997, 1998). 

中 国 大 鲍 (Andrias davidianus) 为 我 国 特有 的 
有 尾 两 栖 动 物 ， 是 世界 上 体形 最 大 的 两 栖 类 物种 
(Song，1982)。 该 物种 在 我 国 所 有 有 尾 两 栖 动 物 中 
分 布 最 广 ， 遍 及 长 江 、 黄 河和 珠江 三 大 水 系 ， 涉 及 
17 个 省 市 的 大 部 分 区 域 (Zhao ，1998; Zhang et al， 
2002; Wang et al，2004)。 但 由 于 多 方面 的 原因 ， 
该 物种 的 栖息 地 质量 急剧 下 降 ， 并 呈现 出 严重 的 片 
断 化 和 破碎 化 (Hu, 1987; Zhang & Wang, 2000; 
Wang et al, 2004). 

分 子 系统 地 理学 是 研究 种 内 地 理 分 布 的 历史 成 
因 、 进 化 历程 和 种 内 遗传 关系 的 重要 途径 。mtDNA 
已 被 广泛 地 用 于 探讨 种 群 差异 与 生物 地 理 的 关系 
(Frankham et al，2002)， 其 中 进化 速度 快 的 D-loop 
(控制 区 ) 是 研究 种 内 遗传 结构 、 遗 传 多 样 性 、 瓶 
颈 效 应 等 最 有 效 的 区 域 (Rosel & Block, 1996; 
David & Tumer，1998)。 在 两 栖 动 物 中 ， 通 过 对 大 
鲍 的 Cyt b 的 比较 发 现 安徽 种 群 具有 独特 的 遗产 特 
ÎE (Muphy et al，2000)， 而 通过 D-loop 基因 探讨 
大 鲍 种 群 遗 传 结构 和 地 理 分 化 ， 其 结果 会 如 何 ? 为 
此 ， 我 们 对 来 自 于 广西 、 河 南 、 陕 西 和 湖南 的 大 鲍 
进行 了 D-loop 基因 全 序列 的 分 析 比 较 研 究 。 


1 材料 和 方法 


1.1 样品 来 源 

本 次 研究 用 的 28 个 大 鲍 样 品 (肌肉 组 织 )， 分 
别 来 自 于 广西 龙 胜 (11 个， 珠江 单 元 )、 河 南 卢 氏 
(11 个 ， 黄 河 单元 )、 陕 西汉 中 (4 个 ) RUSEESSETR 
山 (2 个 ) (长 江 单元 ) (图 1)。 样 品 采 于 2002 一 
2004 年 ， 用 75% 无 水 乙醇 浸泡 并 保存 于 -20 CH 


冰箱 中 。 
1.2 总 DNA 提取 

取 约 50 mg 样品 ， 剪 碎 ， 放 人 1.5 mL 的 Ep- 
pendorf 离心 管内 ， 用 蛋白 酶 K 消化 过 夜 ， 然 后 依 
次 用 饱和 酚 : 氯 仿 : 异 戊 醇 (25:24:1) RISUS: FR 
醇 (24:1) 分 别 抽 提 2 次 ， 再 用 无 水 乙醇 ( - 20 
C) 沉淀 和 TPH LEHR, BOAT, REA 
T 40 uL K TE P., 

1.3 PCR } 3$ mtDNA D-loop 基因 全 序列 

扩 增 引物 参考 Zhang et al (2003) 对 大 鲍 mtD- 
NA 全 序列 的 测序 结果 进行 设计 ， 成 为 与 大 鲍 D- 
loop 全 序列 完全 匹配 的 引物 序列 为 DLP; 5- 
GACTCTTCTTTCACTTTTTCA-3’, DLL: 5'-AGAT- 
GTCTTGGGCTATTGT-3' , 

PCR 反应 在 PTC-200 Peltier Thermal Cycler (MJ 
Research, Inc. USA) 上 进行 。 反 应 总 体积 为 50 
4L。 反 应 程序 为 ，95 CHEP 10 min, 94 CEH 
1 min, 58 CHIE 1 min, 72 CHER 1 min， 并 进行 
35 个 循环 ; 72 CEP 10 min。 利 用 未 加 模板 DNA 
的 反应 液 作 为 空白 对 照 ， 以 检查 是 否 有 污染 存在 。 
PCR 产物 用 1% 琼 脂 糖 凝 胶 电泳 检测 ， 于 - 20 CR 
存 。 

1.4 PCR 产物 纯化 和 测序 

扩 增 得 到 的 PCR 产物 经 琼脂 糖 凝 胶 电 泳 ， 所 
需 带 切 下 后 用 宝生 物 工程 (KE) ERAF 
(Takara) 生产 的 DNA 凝 胶 回 收 试剂 盒 (TaKaRa A- 
garose Gel DNA Purification Kit Ver. 2.0) 纯化 回 
收 。 最 后 根据 琼脂 糖 电 泳 检测 和 定量 。 全 部 序列 均 
由 上 海 晶 泰 生物 技术 公司 用 Licor 4300 基因 分 析 仪 
进行 测序 。 

1.5 DNA 序列 数据 的 处 理 

采用 Clustal X 软件 将 测序 结果 进行 对 位 排序 ， 
确定 变异 位 点 及 单 倍 型 。 用 DnaSP4.0 软件 计算 各 
地 理 种 群 的 核 苷 酸 多 样 性 (nucleotide diversity, 
x)、 单 倍 型 多 样 性 (haplotype diversity, h) 及 种 群 
间 的 核 苷 酸 歧 异 度 (nucleotide divergence, Dxy), 
并 进行 Fu 的 Fs 中 性 检验 。 用 Arlequin1.1 中 的 
AMDVA 分 析 方 法 计算 地 理 单 元 间 的 分 化 水 平 CF- 
statistics, Fst) 及 其 基因 流程 度 (Nm )， 以 进 一 
步 揭示 群体 的 分 化 程度 .用 MEGA2.0 软 件 中 的 Kim- 
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图 1 本 文大 鲍 样 品 的 采样 地 (黑色 圆 点 表示 ) 和 Murphy 的 采样 地 〈 正 方形 表示 ) 
Fig. 1 Sampling sites of the Chinese giant salamander in this study (the dots marked 
in black) and Murphy's study (the rectangle marked in black) 
闭合 线圈 表示 大 鲍 目 前 的 分 布 范围 (Wang et al, 2004). 


Cycles are the recent distribution areas of Chinese giant salamander (Wang et al, 2004). 


ura (1980) 双 参 数 法 (2-parameter) 分 析 了 各 单 倍 
型 之 间 的 遗传 距离 。 用 MEGA2.0 软件 中 的 邻近 距 
离 法 和 PAUP4.01 最 大 简约 法 构建 系统 进化 树 ， 均 
以 新 疆 北 鲍 ( Ranodon sibiricus) 作为 外 群 ， 同时 
应 用 自 举 检验 (bootstrap test) 估计 系统 树 中 节点 
的 置信 度 。 


2 结果 与 分 析 


2.1 mtDNA D-loop 的 全 序列 及 其 差异 

在 28 RAKE mDNA D-loop 771 bp 的 全 序列 
H, HEM 64 个 多 态 性 核 萌 酸 位 点 ， 占 全 部 碱 基 
数 的 8.2690; 简约 信息 位 点 19 个 其 中 转换 
(transition) 6 个 和 其 换 (transversion) 2 个 , TA 
缺失 (insertion/deletion) 11 个 。 序 列 分 析 结 果 表 
Hj, 在 4 个 地 理 种 群 的 28 个 样品 中 共存 在 单 倍 型 
27 种 ; 广西 种 群 10 种 ， 河 南 种 群 11 种 ， 湖 南 2 
种 ， 陕 西 4 种 。 同 时 根据 这 4 个 地 理 种 群 的 地 理 分 
Tg (图 1)， 分 成 3 个 地 理 单元 ， 即 珠江 单元 ( 广 
西 种 群 )、 黄 河 单元 (河南 种 群 ) 和 长 江 单元 (W 
南 和 陕西 种 群 )。 
2.2 ”地理 单 元 内 和 单元 间 的 遗传 结构 

在 27 个 单 倍 型 中 ， 各 地 理 单元 之 间 并 没有 共 


享 的 单 倍 型 ， 且 单 倍 型 间 的 序列 差异 在 0.1% ~ 
2.7% ， 平 均 为 1.32% 。 珠 江 单元 、 黄 河 单元 和 长 
江 单 元 的 单 倍 型 多 样 性 指数 (A). 分 别 为 0.00150、 
0.00215 和 0.00226， 而 核 苷 酸 多 样 性 指数 (x) 为 
0.00913、0.01114 和 0.00936。3 个 单元 的 单 倍 型 
多 样 性 指数 和 核 苷 酸 多 桩 性 指数 的 值 都 偏 低 ， 而 且 
珠江 单元 的 这 两 个 指数 值 低 于 长 江 和 黄河 两 个 单元 
WÉ CE D. M Fs 检验 结果 (Fs = - 17.054) 表 
明 ， 大 鲍 种 群 经 历 了 一 个 强 的 种 群 增长 和 扩 撤 ， 近 
期 存在 遗传 变异 [根据 Fu(1997)、Su et al(2001) 
Fs 趋向 于 负 值 的 推断 ] 。 

根据 分 子 变异 分 析 (AMOVA) 结果 显示 C 
2) ， 珠 江 单元 和 黄河 、 长 江 单元 之 间 分 化 程度 显著 
(P<0.001)， 而 长 江 单元 和 黄河 单元 之 间 差 异 不 
显著 (已 > 0.05)。 同 时 地 理 单元 内 分 化 程度 占 
99.31%， 而 单元 间 只 有 0.69% ， 表 明 遗 传 差异 主 
要 发 生 在 单元 内 ， 而 且 各 地 理 单元 之 间 的 基因 流 较 
频繁 。 因 此 ， 珠 江 单元 的 广西 种 群 与 其 他 地 理 种 群 
有 显著 的 遗传 差异 。 
2.3 系统 地 理 分 布 格局 关系 

构建 的 NJ 树 和 MP 树 (图 2) 显示 ，27 个 单 
倍 型 散布 在 不 同 的 分 布 群 中 ， 呈 现 一 种 混杂 的 分 布 
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格局 ， 并 未 分 成 代表 3 个 地 理 单元 (长 江 单元 、 黄 
河 单元 和 珠江 单元 ) 的 聚合 群 。 
3 讨论 

一 般 认 为 ， 两 栖 类 动物 的 D-loop 基因 长 度 普 
遍 长 于 哺乳 类 和 鸟 类 (Roe et al, 1985; Zardoya & 
Axel, 2000), EKRA 771 bp， 明 显 小 于 其 他 
两 栖 类 、 哺 乳 类 和 鸟 类 。Wenink et al (1994) iA 
为 ， 两 栖 类 D-loop 是 整个 mtDNA 分 子 中 进化 最 快 
的 一 个 区 域 ， 比 mtDNA 中 的 其 他 区 域 快 3 ~ 5 倍 。 
然而 ， 根 据 我 们 对 大 鲍 D-loop 基因 全 序列 的 差异 
分 析 发 现 ， 在 大 鲍 的 27 个 单 倍 型 间 ， 其 平均 变异 
只 有 1.32%， 稍 高 于 Cy b 基因 序列 的 变异 
(1.290) ( 待 发 表 ) ， 同 样 有 别 于 两 栖 类 的 其 他 物种 
(Yang et al，1994)。 也 说 明了 大 够 的 遗传 多 样 性 
比较 贫乏 。 因 此 ， 大 鲍 是 两 栖 类 中 比较 独特 的 物 
种 。 

已 有 研究 (Hu, 1987; Zhang et al, 2002; 
Wang et al, 2004) H, ASTRI E Hi H EDI RUE E 
A3E B] S6, URK A K EEE T 
片断 化 的 格局 ， 呈 现 出 大 小 不 等 的 11 个 分 布 板块 
(图 1)。 我 们 大 鲍 样 品 的 湖南 、 陕 西 、 河 南 采 集 地 


和 Murphy et al (2000) 湖南 、 陕 西 、 山 西 、 河 南 
与 四 川 的 样品 采集 地 都 同属 于 一 个 大 分 布 区 (图 
1) ， 而 广西 和 安徽 样品 采集 地 却 分 别 属 于 两 个 单独 
的 小 区 域 。Barber (1999) 认为 地 理 障 碍 在 种 系 地 
理学 结构 的 形成 过 程 中 可 能 起 到 非常 重要 的 作用 。 
在 空间 跨度 上 ， 广 西 种 群 离 中 心 的 大 块 分 布 区 较 
近 ， 而 安徽 种 群 较 远 ,但 都 存在 地 理 隔离 。 

Feller & Blair (1998) 认为 两 栖 类 物种 亲缘 关 
系 的 形成 与 历史 地 理 变化 有 极 大 的 关系 。 从 历史 地 
理学 上 看 ， 距 今 约 150 百 万 年 的 中 生 代 侏 罗 纪 ， 印 
支 运 动 的 作用 力 形成 了 原始 的 云贵 高 原 ; 约 100 A 
万 年 的 白垩 纪 发 生 的 燕山 运动 使 云贵 高 原 上 升 ， 成 
了 长 江 和 珠江 的 分 水 岭 ， 并 使 珠江 流域 产生 一 系列 
断裂 破碎 带 ， 以 后 形成 的 珠江 水 系 蚁 昨 于 这 些 断裂 
带 中 (Situ，2001)。 在 40~ 30 百 万 年 前 的 新 生 代 ， 
喜马拉雅 运动 使 青藏 高 原 出 水 成 陆 ; 之 后 形成 了 系 
统 的 长 江水 系 (Lu & Chen，2001)， 而 约 3 百 万 年 
前 的 第 四 纪 是 黄河 形成 、 发 育 和 演变 的 地 质 时 代 
(Chen，2002)。 可 见 ， 珠 江 的 形成 早 于 长 江 和 黄 
河 ， 又 由 于 有 云贵 高 原 的 阻碍 而 与 长 江 和 黄河 分 
离 。 所 以 我 们 认为 由 此 而 产生 的 早期 大 鲍 种 群 应 分 
化 成 珠江 单元 和 长 江 、 黄 河 单元 ,并 独自 演变 形成 


R1 大 鳃 3 个 地 理 单元 的 核 背 酸 歧 异 度 (Dxy， 对 角 线 上 , %)、 单 倍 型 多 样 性 (5) 和 核 苷 酸 多 样 性 (m) 
Tab. 1 Nucleotide divergence (Dxy, above diagonal, %), haplotype diversity (h) and nucleotide diversity (1) in three 


geographic units of Chinese giant salamander 
珠江 单元 黄河 单元 


Pearl river region Yellow river region 


Yangtze river region 


珠江 单元 
Pearl river region 


黄河 单元 


Yellow river region 
长 江 单元 


Yangtze river region 


1.050 + 0.008 


长 江 单元 单 倍 型 多 样 性 (h) AEREE (x) 
Haplotype diversity Nucleotide diversity 
0.963 + 0.006 0.00150 0.00913 
1.027 + 0.007 0.00215 0.01114 
0.00226 0.00936 


表 2 X8 3 个 地 理 单元 之 间 分 化 程度 (Fu, HARF) HAH (Nm, HAREL) 以 及 遗传 变异 在 地 理 单元 中 的 分 布 
Tab. 2 F, (below the diagonal) and Nm (above the diagonal) among the three geographic units of Chinese giant sala- 
mander, and distributions of genetic variation among the regions and genetic variation within the regions 




















珠江 单元 黄河 单元 长 江 单元 
Pearl river region Yellow river region Yangtze river region 
珠江 单元 Pearl river region 54.505 无 
黄河 单元 Yellow river region 0.00909 *** 50.261 
长 江 单 元 Yangtze river region o" 0.00985 
地 理 单元 间 的 遗传 变异 Genetic variation among the regions. (96) 0.69 
地 理 单元 内 的 遗传 变异 Genetic variation within the region (96) 99.31 


“”“”P<0.001 (AMDVA). 
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HEN? 
GX8 
GX5 SX4 
HEN6 HEN6 
GX9 
Hs HEN8 
HEN! GX1 
HENS 66 HUN2 
GX! Sxi 
66 HUN2 HUN! 
92 SX4 79 GX9 
HUN! GX5 
50 GX10 HENIO 
sx! GX10 
GX2 GX2 
HEN2 HEN9 
72 HEN9 68 HEN2 
HENIO GX4 
SX3 GX3 
53 CI 67 SX3 
cx: GX7 
89 HEN7 96 HEN? 
HENS HENS 
HEN3 HEN3 
7l GX6 37 GX6 
76 SX2 80 SX2 
HEN4 HEN4 
Ranodon Ranodon 
ANJ 树 B:MP 树 


图 2 用 NJ 法 (A) 和 MP 法 (B) 构建 的 大 够 28 个 地 理 种 群 mDNA D-loop 基因 序列 的 系统 树 
Fig. 2 Phylogenetic tree among 28 populations of Chinese giant salamander using NJ method (A) and 
MP method (B) inferred from mtDNA D-loop gene sequences 
GX: 广西 种 群 (Cuangx population); SX: 陕西 种 群 (Shaanxi population); HEN: 河南 种 群 (Henan population); HUN: W 


南 种 群 (Hunan population); Ranodon: 外 群 (Outgroup )。 


一 定 的 遗传 结构 差异 。 此 外 ，Jiang (1998) 认为 ， 
约 5 百 万 年 前 的 晚 第 三 纪 以 来 发 生 新 结构 运动 的 断 
裂 抬升 ， 造 成 了 安徽 黄山 海拔 较 高 的 地 势 ; 从 而 形 
成 大 鲍 安 徽 种 群 的 一 个 独立 山区 ， 更 易 发 生 遗 传 漂 
AE, 导致 该 种 群 与 其 他 地 理 种 群 的 遗传 分 化 《Mur- 
phy et al, 2000), 

但 是 ， 本 文 构建 的 NJ 树 和 MP 树 并 未 显示 出 
推定 的 3 个 地 理 单元 ， 且 各 种 群 之 间 的 基因 流 也 比 
较 频 繁 。 推 测 这 可 能 是 现 有 的 地 理 隅 离 格 局 的 形成 
时 间 还 不 够 长 ， 各 种 群 之 间 还 未 达到 具有 显著 遗传 
隔离 和 分 化 的 单 系 水 平 ; 或 是 近期 非常 严重 的 人 为 
干扰 所 致 。 因 此 ， 鉴 于 目前 大 鲍 的 破碎 化 分 布 ， 建 
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